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We examine for standard-like orbifold compactification models the constraints due to 
quarks and leptons generation non-universality of soft supersymmetry breaking interac- 
tions. We follow the approach initiated by Ibanez and Lust and developed by Brignole, 
Ibanez and Munoz . This is based on a locally supersymmetric a— model action of 
moduli and matter fields obeying the stringy duality symmetries. It is assumed that the 
low energy fields of the minimal supersymmetric standard model are in one-to-one cor- 
respondence with string massless modes and that supersymmetry breaking takes place 
simultaneously with the lifting of flat directions for dilaton and compactification moduli 
fields. The breaking of supersymmetry is represented in terms of dilaton and moduli 
auxiliary field components and, consistently with a vanishing cosmological constant, is 
parametrized in terms of the dilaton-moduli mixing angle 9 and the gravitino mass scale 
m g . The soft supersymmetry breaking interactions (gaugino masses, squarks and slep- 
tons mass matrices, scalars interactions A and B coupling constants) are calculable as a 
function of these parameters and of the discrete set of modular weight parameters spec- 
ifying the modular transformation properties of the low energy fields. To establish the 
flavor changing neutral current constraints we solve the renormalization group one-loop 
equations for the full set of gauge, Yukawa and supersymmetry breaking coupling con- 
stants. A simplified version is used in which one treats the contributions from the third 
generation quarks and lepton Yukawa couplings exactly, while retaining for the first and 
second generations couplings only the leading order term in the large logarithm variable. 

Numerical results are obtained for the quantities Amn = ^m^mn^ni corresponding to 
the mass matrices of squarks and sleptons in the super-CKM basis, for which experimental 
bounds can be determined via the super-box and super-penguin diagrams with gluino or 
neutralino exchange. 
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1. Introduction 

The string-based 4-d supergravity cr-models inherit from the duality symmetries of 
compactified strings (see ref.[l] for a review) an additional unbroken symmetry corre- 
sponding to discrete modular groups which act as Kahler transformations on the moduli 
and matter superfields [2] . The simplest case where a general analysis of the moduli space 
structure and of the modular symmetry action can be made is that of compactifications 
described by world sheet cr-models with (2,2) supersymmetry [3]. No results of compara- 
ble rigor exist for the phenomenologically more interesting (0,2) models, with or without 
Wilson lines. However, extensive analyses of these models exist within the orbifold toroidal 
compactification framework [4-8]. Here, string tree-level correlators are used to identify 
the Kahler potentials of the string massless moduli and matter modes and to determine 
modular transformation laws of the associated fields leading to assignment of modular 
weights [8]. String one-loop correlators [9,10] are also required to establish the cancella- 
tion of a— model duality anomalies which takes place as a result of string mass threshold 
effects on gauge coupling constants in combination with an analog of the Green-Schwarz 
mechanism [11,12]. 

The restrictions placed by duality on supergravity a-models greatly enhance the pre- 
dictive power of studies of supersymmetry breaking within the so-called standard-like 
models, which are based on the assumption of a direct compactification to the minimal 
supersymmetric model. One important aspect is to narrow the choice of the functional 
dependence on the compactification moduli of the non-perturbative component of the su- 
perpotential which is supposed to lift flat directions in dilaton and moduli field space. The 
main hope, however, is that these non-vanishing VEVs (vacuum expectation values) of 
scalar fields components of dilaton and moduli are accompanied by non-vanishing VEVs 
for the corresponding auxiliary fields components, thus associating supersymmetry break- 
ing with the same (as yet unknown) mechanisms that fix the gauge and Yukawa coupling 
constants of the low energy theory [13-15]. The viability of this picture has been exam- 
ined in several recent works within the hidden sector gaugino condensation approach with 
a view to make contact with phenomenology [16,17]. However, to avoid occurrence of a 
vanishing dilaton auxiliary field component, Fs = 0, or a non- vanishing cosmo logical con- 
stant, elaborate versions must be used, in which hidden matter condensation or gaugino 
condensation from several gauge group factors play a role in the local supersymmetry 
breaking [18]. Motivated by these limitations, Brignole, Ibahez and Muhoz [19] have pro- 
posed an alternative model independent approach avoiding any specific assumption on the 
non-perturbative superpotential except for the condition of non-vanishing auxiliary fields 
components of dilaton, moduli and possibly matter superfields, which are parametrized in 
such a way as to ensure a vanishing cosmo logical constant. While still rooted within the 
gaugino condensation framework, this approach is characterized by a stronger emphasis 
on auxiliary fields over scalar fields, thus providing a transparent representation of super- 
symmetry breaking in terms of mixing angles within the goldstino and the gravitino mass. 
However, the scalar VEVs stand out as extra free parameters in this approach. 

The low energy limit of standard-like models can then be derived on minimal grounds 
by assuming a one-to-one correspondence between string massless modes and the low en- 
ergy (quarks, leptons, Higgs) particles and, in the familiar way, expanding the supergravity 
a— model model action in powers of the inverse Planck mass 1/Mp. The supersymmetry 
breaking parameters (gaugino masses, scalars masses, bilinear and trilinear scalars inter- 
actions couplings) can be explicitly calculated in terms of the dilaton and moduli scalar 
and auxiliary components VEVs and of the fields modular weights. 

One characteristic feature of the predicted soft supersymmetry breaking coupling 
constants is the suggested presence of large non-universal contributions with respect to 
the quarks and leptons generation quantum numbers [81. Such generation dependence of 
squarks and sleptons mass matrices and of A-matrices for couplings of squarks and slep- 



tons to Higgs bosons could contribute to flavor changing neutral current processes at an 
observable level, a fact which can be used under appropriate assumptions to set upper 
bouds on generation dependence [20,21]. However, the relevant physical observables in 
supergravity models are those obtained after summing the large logarithms contributed 
by radiative corrections. This requires solving the renormalization group flow equations 
down to the electroweak scale, using the soft parameters as boundary conditions at the 
string unification scale [22,23]. While the generation independent interactions of gauginos 
reduce whichever large non-universality present at the string scale, Yukawa interactions 
have an opposite effect. It is then mainly a quantitative issue to determine what residual 
non-universal dependence is not wiped out at the electroweak scale. Also, the fact that 
the quantities which enter in the evaluation of super-box or super-penguin diagrams in- 
volve the tranformation matrices which diagonalize the quarks and leptons mass matrices, 
implies that the relevant observables are also sensitive to the flavor structure of Yukawa 
interactions. 

The issue of universality of supersymmetry breaking is closely tied with the fla- 
vor problem. Conceptually, it is difficult to reconcile the important flavor asymmetries 
present in Yukawa interactions of quarks and leptons with the postulate of a spontaneous 
supersymmetry breaking independent of flavor or generation quantum numbers [24] . Non- 
universality also poses an acute naturalness problem to supergravity models [25,26]. This 
arises in grand unified theories embedding of the minimal supersymmetric standard model 
where radiative corrections from decoupling of superheavy modes can induce contribu- 
tions to mass matrices of scalar superpartners in generation space which are misaligned 
with those arising from Yukawa couplings of low energy modes [20]. The early envisaged 
possibility to limit the size of non-universal terms on the basis of flavor (horizontal) dis- 
crete or continuous symmetries has been revived recently in a number of interesting works 
[27-30]. Relaxing the universality restriction also plays an important role in widening the 
allowed parameter space of the minimal supersymmetric standard model [311. Note how- 
ever that string theory interpretation of non-universality differs from that or field theory 
in two respects: (i) Non-universality is associated with generation dependent assignment 
of modular weights to low energy fields taking discrete rational values rather than with 
the spontaneous breakdown of an horizontal symmetry; (ii) It sets in at the larger string 
unification scale which should ease its smearing by flavor-blind radiative corrections. Thus, 
string theory provides us with a framework to accommodate large non-universality as well 
as mechanisms to reduce their observable effects. Related recent works on this subject are 
refs. [32,33] 

Our main goal in this work is to examine the implications of the experimental fla- 
vor changing neutral current bounds [20,34] concerning the presence of a generation non- 
universality in the dimension-two scalars mass matrices, rh 2 MNl [M,N = L,R], and the 
dimension-three scalars couplings matrices Afj, [x = u, d, e\. In the super-CKM basis, with 

flavor diagonal Yukawa couplings of quarks and squarks or leptons and sleptons to gaugi- 
nos, the generation dependence in the super-box or super-penguin diagrams, propagating 

superpartners of low energy particles, involve the quantities, Amn = Vm^'mm^n-* where 
Vm, [M = L, R] are transformations which diagonalize the quarks and leptons mass matri- 
ces. We shall restrict consideration to the standard-like orbifold compactification leading 
to the minimal supersymmetric standard model . The local supersymmetry breaking pa- 
rameters are then determined as a function of the dilaton and moduli auxiliary and scalar 
fields components at the large string unification scale. To compare with upper bounds on 
the parameters Amtv, we need to integrate the scale evolution equations. As is well known, 
the corrections from gauge bosons and fermions loops, being flavor-blind, cause important 
dilution of flavor mixing effects [26]. Thus, non-universality should be strongly suppressed 
in a dilatonic goldstino vacuum, cos# ~ [19]. A naive estimate [19] including only gauge 
interactions in the scale evolution yields a large upper bound on the cosine of the mixing 



angle, | cos#| < 1/V3, corresponding to a narrow range of values for the dilaton- moduli 
mixing angle, 90° < 6 < 55°. An independent source of flavor mixing is also present in the 
trilinear scalar interactions, which are affected too by the renormalization group evolution. 

The contents of the paper are organized in 4 sections. In Section 2, we recall the main 
ingredients of supergravity a— model effective actions incorporating duality symmetries, 
as developed by Ibanez and Lust [8] ( Subsection 2.1), and of supersymmetry breaking in 
this framework, using the parametrization of Brignole, Ibanez and Muhoz [19] (Subsection 
2.2). Next, we present the approximate version used for the one loop renormalization 
group equations in which the contributions from the third generation particles (top and 
bottom quarks and r-lepton) are included exactly, while those of the first and second 
generations are included perturbatively (Subsection 2.3). Finally we discuss the choice at 
string scale for the Yukawa coupling matrices (Subsection 2.4). In Section 3, we present 
and discuss our numerical results. In Section 4, we summarize our main conclusions and 
suggestions for a future improved treatment of the problem. 

2. Duality symmetry and supersymmetry breaking interactions 

2.1 a— model effective action 

By moduli space(s) of a superstring theory one means the configuration space sector(s) 
of the associated world sheet field theory spanned by marginal operators deformations 
preserving the superconformal symmetry. This space constitutes the vacuum manifold 
of the low energy 4-d effective field theory. The moduli space of string-based <7-models 
supergravity are parametrized by: (1) the space-time dilaton field, S(x); (2) the compact- 

ification moduli fields, M(x), M(x), associated to deformations of the Kahler and complex 
structures of the 6-d internal space manifold; (3) the matter-like background fields, A^ , as- 
sociated to Wilson gauge flux lines winding around non-contractible loops in the lattices of 
Eg x Eg root space and internal 6-d space-time and other possible massless fields associated 
with duality [6]. Orbifolds (see refs.[35] for reviews) O = R 6 /(AxP) represent special exact 
solutions for an internal space given by a 6-torus of lattice A with a discrete point symmetry 
group P generated by elements of the SU(3) holonomy subgroup of the 6-d tangent space 
group SO (6). The point groups of abelian, Zn or Zn x Zm, orbifolds are generated by ele- 
ments which can be represented in a complex basis as: 9 = diag{e 2lT%di )^ [£V 9i = 1] with 
complex conjugate indices i, i — 1, 2, 3 labeling coordinates, X l=1,1 = -4=A M=1±l2 , • • •, and 

spinors. The relevant compactification moduli fields constitute then a subset of the mod- 
uli matrices, [T^, Tfj], [U^, U{j], [i,j = 1, 2, 3], invariant under P. For convenience, we 
shall restrict ourselves to the so-called generic orbifold case, characterized as allowing only 
diagonal Kahler structure moduli, Ti,Ti[i = 1,2,3], or a subset of these, while excluding 
entirely complex structure moduli Ui, Ui and matter background fields A^ . (See ref.[9] for 

the list of allowed moduli and ref.[8] for the list of a— model moduli spaces). As in most 
of the recent studies, we also shall specialize to decomposable 6-d tori with a metric tensor 
given by a direct product of metric tensors for three orthogonal 2-d tori, or complex planes, 
labelled by indices i, i = 1, 2, 3. (Calculations with non-decomposable tori are considered 
in refs. [36,37]). 

The Kahler potential of the effective N = 1 supergravity effective action can be 
derived in general form for the case of world sheet with (2,2) supersymmetry [3]. The 
part involving untwisted moduli is given in exact form as: 

K moduli (S, S, M, M, A™, <t) = _ J2 log U + Tl - J2 < f II e^ v ^) A A 

4 ^ -\i n ' 



-^log(^ + ^)-log(5 + 5), (1) 

i 

while the part involving untwisted and twisted sectors matter superfields, denoted generi- 
cally by A a = A a + Q\ a + 6 2 F a , where the index a designates whatever quantum numbers 
(gauge, family,...) are needed to characterize string modes, is given in an expansion in 

powers of jf- as: 

K matter (M, M,A a , A^) = ]jT nm +^ ^ 

a i m ^ a ' 

(2) 
where n l a , I™ are fields modular weights, rational numbers characterizing the modular group 
transformation laws. For each matter field, we have inserted its minimal gauge coupling 
interaction, such that each group factor G a , with coupling constant g a , is associated the 

real superfield V a (R a ) = —iV£ T a (R a ), with T a (R a ), the Lie algebra generators in the 
representation R a subject to the normalization convention Tr(T a T a i) = \c{R a )5 aa ' , where 
c(R a ) is the Dynkin index of representation R a . Note that the absence of off-diagonal 
kinetic terms in eq.(2) is an exact result for (2,2) models which is expected to remain true in 
general cases. The above formulas are found to hold for the wider class of world sheet (0,2) 
models only for orbifolds. We shall restrict ourselves to this case in the sequel. The duality 
symmetries are realized for the diagonal moduli, T^ = Ti, by the modular symmetry group 

T = [SL(2, Z)\ (or an appropriate subgroup), whose action on the moduli and matter 

fields in the supergravity basis (as opposed to the string vertex operators basis ) [6] is 
specified by the transformations: 

A a ^Y[(ic l T l + d l )<A ct , (3) 

i 

= 1, 2, 3, define elements of SL(2, Z)t { and 
the modular weights are given as [8]: n l a = —5 a i, (untwisted sector a = i = [1,2,3]) and 
n l a = — (1 — 9 l + p l — (f ), (twisted sectors), p\ q l being the number of complex conjugate 
left-moving sector oscillator excitations of type cx l rn _ ei , a % m+ei . Note that the background 
fields A^ transform as untwisted states and that for complex structure moduli analogous 
formulas hold for the associated modular weights I™ with the interchange of p 1 with q l . 

Each gauge symmetry group factor G a of the full gauge group, ria^«' ^ s associated 
a holomorphic gauge function f aa >, defined in the adjoint representation. In the glob- 
ally supersymmetry limit, with the field strength chiral superfield defined as, —iW£T a = 

jD 2 e~ Va D ce e Va = —\{i\% + • ■ ■)T a 1 this contributes the familiar Lagrangian: 

Lgauge = \j^ [ ^^^ W * f**' & TJ + C.C., (4) 

such that f aa '(S, Mi) = k a S5 aa ' , with the dilaton VEV being related to the string coupling 
constant as < S >= l/g\ and k a being the level of the world sheet Kac-Moody algebra 
for the group G a . The only requirement from duality on the superpotential is that it be a 
modular form of SL(2, Z)^ of weight —1, corresponding to the modular group acting as the 

Kahler transformation, K — >■ K + fi + fi,W — *■ e~^W, [fi = \og(iciTi + d{), fa = /*]. 
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The consideration of string and a— model loop expansions provides important in- 
formation on the moduli dependence of the effective action. Realization of duality as a 
Kahler transformation, entails the presence of chiral 11(1)^ R-symmetries transforming 

matter and gauge fermions as, ip — > ipe~ tWi ^' Im ^ i , with weights Wi(X a ) = |, Wi(xa) = 

— 1(1 + 2n l a ). For the diagonal modular subgroup, / = ^2 i fi,w(X a ) = \,w(xa) = 
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At the quantum level these symmetries are affected by 



a— model triangle anomalies involving non-local couplings of gauge bosons to Kahler and 
matter composite connections, 



jYl [ d 2 0W a W a 
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(Analogous a— model gravitational anomalies are also present.) The duality anomalies 
contain a gauge group independent (universal) part which can be cancelled by an ana- 
log of the Green-Schwarz mechanism [11,12]. This is implemented in the dilaton chi- 
ral multiplet formulation by changing the dilaton field variable to [8,12]: S — > S new = 
S — ,-Aa Yli $hs l°g(^i + Ti), such that the new field transforms under modular transfor- 
mations as: S new — ► S new + ,\i Y^i $Gs(fi ~~ /*)> wnue $ is inert. Performing the implied 

substitution in eq(l), and suppressing the subscript 'new', the corrected dilaton Kahler 
potential becomes: 



log(5 + S) -> log(y) = log (s + S + j^ J2 6 gs Iog(T 4 + ?A . 



(5) 



where the modular invariant function Y is interpreted as an effective dilaton field whose 
VEV defines a renormalized string coupling constant field, < Y >= 2/g\. At string one- 
loop level, there also arises mass threshold corrections in the gauge coupling constants 
which, after changing the dilaton field variable in the tree level term, leads to corrected 
gauge functions [9]: 



faa'{S,Ti) = f a (S,Ti)5 aa > = 5 aa > 



w+E^^cr*) 



(6) 



where rj(T) is Dedekind function (automorphic modular form of weight |) and 6 a l = 

ba ~ k a S GS , with b^ = ^[c(G a ) — ^2 a (l + 2n l a )c(R a )], corresponding to the /^-function slope 
parameters for the N=2 suborbifolds associated with the subgroup of the point symmetry 
group leaving the i— plane fixed. Using eqs.(l) and (2) to write the above non-local effective 
Lagrangian for the a— model modular anomalies as: 

l2n2 



~Ej A "" r ^ £ tf MT, + T t ) + c.c, 



we see that the modular transform of L n i is composed of a moduli-dependent term which 
is cancelled by threshold corrections (second term in the gauge function, eq.(6)) and a 
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universal moduli-independent term which is cancelled by virtue of the Green-Schwarz 
dilaton- moduli mixing (first tree level term in eq.(6)). Parenthetically, we note that an 
analogous cancellation also takes place for the QCD triangle anomalies affecting the Kahler 
R-symmetries 11(1)^ which are compensated by the breaking term occurring in threshold 

corrections, {9 a ) thres = -(4n) 2 Imf a (S, Tj) = -b'* log ( |^ J . This term arises [37] be- 
cause any modular invariant coupling term quadratic in colored fields can be transfered 
via a chiral rotation into a contribution to the ^-vacuum parameter, whose variation under 

U(l) Ti , Si9 a = Uc(G a ) - E Q (1 + 2n l a )c(R a ) )(fi ~ fi), is exactly canceled by that of 

threshold corrections, 5i(9 a ) thres = -b'*(fi - fi). 

To simplify the discussion of moduli dependence, it is convenient to limit the number 
of unknowns by restricting to a one-dimensional direction in moduli space, corresponding 
to the overall volume modulus field T(x). For an isotropic orbifold with three equivalent 2- 
tori, the overall modulus is defined by the restriction: T\{x) = T 2 (x) = T^(x) = T(x), and 

transforms under the diagonal subgroup SL(2, Z) of [SL(2, Z)\ , whose modular weights 
are given by the sums, n a = J^ i n l a . The overall modulus for an anisotropic orbifold can 
be defined analogously by the condition: r{T\ = r 2 T 2 = r^T^ = T, where r$ are real 
parameters characterizing the geometrical shape of the 6d-torus which obey Y\.i r i = 1? by 
virtue of volume conservation. The squeezed orbifold with T\ ^> T 2 , T3 has r± ^ir 2l r^. The 
diagonal moduli VEVs are related to the radii of the 2-d tori as: Re(Ti) = cnR 2 /((47r) 2 )o:'), 
with a' the string slope parameter and c 2 N oc det(G), with G the 6-torus metric tensor. The 
cn are calculable constants depending on the Zn orbifold order. For instance, C3 = cq = 
\/3, C4 = 2\/2. (Applications dealing with non-decomposable 6-d tori are developed in 
refs.[36].) Naturally, the values at which these radii, or the corresponding overall modulus 
and anisotropies, settle is determined by compactification dynamics which lies outside 
the present framework. However, because compactification means a choice of vacuum, 
minimization of the supergravity scalar potential will place certain conditions on these 
VEVs. Thus, for decomposable tori, the equivalence of different planes imposes a discrete 
permutational symmetry among the allowed moduli. An anisotropic solution could then 
arise in either of the following two situations: (1) A subset of the diagonal moduli is allowed. 
The Za, orbifold, for example, allows for two diagonal moduli fields only; (2) Spontaneous 
breakdown of the permutational symmetry. Nevertheless, according to ref.[14], it is not 
easy to find generic moduli configurations which realize the latter possibility. 

Let us briefly summarize the main conclusions of Ibanez and Lust [8] concerning the 
range of values that can be assigned to modular weights. Given an orbifold, with point 
group in the allowed list of abelian Zn or Zn x Zm groups [38], this is associated dis- 
crete sets of twists 0\ , [p = 1, • • • , T] oscillators modings, Kac- Moody levels, and fixed 

points f\ ; , [q = 1, • • • , F]. While Yukawa couplings depend on all these items, modular 
weights are essentially fixed by the first three only. The larger are the values assigned 
to k a and to oscillator excitations, the wider is the range of modular weights. The main 
issue for minimal supersymmetric standard model compactification is to assign modular 
weights to the low energy fields consistently with cancellation of duality anomalies and 
with gauge coupling constants unification. An important characteristic of orbifolds here 
is the number of overall unrotated planes (i.e., planes left fixed by sweeping through the 
entire set of twists). This number can take the values: N unr = 3,2, 1,0. The orbifolds 
Z N x Z M have N unr = 3; ZL has N unr = 2; Z 4 , Z 6 , Z 8 ,Z' 8 , Z 12 , Z[ 2 all have N unr = 1; 
prime number orbifolds, Z3, Z7, have N unr = 0. The duality anomalies cancellation places 



the constraint, 6 a * = 0, for overall rotated planes, i = 1,2,3, where the index a labels 
the standard model group factors, a = 3,2, 1, corresponding to SU(3) x SU(2) x U(l). 
For the minimal supersymmetric standard model the conditions b^ = for an overall 
rotated plane, take the form: 6 a l = — c a — J2 q c^n l a = k a 5 GS , where g is a genera- 
tion index, cs = 3, c 2 = 5,ci = 11 and c" are rational numbers, which for standard 
model fields labeled in the following order, a = [Q, U c , _D C , L, E c , Hi, H 2 ], are given by 
[8]: cf = [2,1,1,0,0,0,0], c£ = [3,0,0,1,0,1,1], c? = [§,§,§,1,2,1,1]. These mod- 
ular anomalies cancellation conditions are absent for the Z^ x Zm orbifolds, because of 
N U nr = 3, and involve one equation of the above type (for each a = 3, 2, 1) for the Z' 6 
orbifold, two for all other non-prime Zn orbifolds and three for prime orbifolds. 

The 5 GS are extra parameters which are calculable on a model-by-model basis, for 

example, by identifying 6 a * = 6 a * — k a d GS to the slope parameter of the N = 2 suborbifold 
generated by subgroups of the point symmetry group leaving the i-plane fixed. Since 
threshold corrections are non-vanishing for overall unrotated planes only, it follows that 
the unification constraints act in complementary manner to the anomalies cancellation 
constraints, yielding stronger conditions for orbifolds with larger N unr . Of course, it is 
always possible to partially relax these constraints by considering anisotropic orbifolds, 
since a condition such as T\ ^> T 2 , T3, which freezes T 2 , T3, renders ineffective any threshold 
corrections involving these moduli. The analysis of Ibanez and Lust [8] indicates that most 
Zjy orbifolds fail to satisfy the combined anomalies cancellation and gauge unification 
constraints, unless one considers higher Kac-Moody levels for the non-abelian groups, 
which is not a favored option in orbifolds constructions. The acceptable canditates, within 
the minimal framework, are: ZL, Z 2 x Z% and Z 2 x Z' 6 orbifolds. Generalization to squeezed 
orbifolds extends the acceptable orbifolds to the whole Zn x Zm set. Recall here that 
anisotropic compactification could also represent a viable mechanism to generate large 
hierarchies in the Yukawa couplings matrices [39]. Unlike the superpotential and gauge 
functions, the Kahler potential is negligibly affected by anisotropies. Indeed, the constraint 
J2i l°g r i — removes all dependence here on r^ at tree level. At one- loop level, a residual 
contribution may arise from the term, £V 5 l GS log(Tj +Tj) = — ^\ 5 GS logr^ +5gs log(T + 
T), which can be absorbed inside the dilaton field. However, whereas threshold effects 
involve components 5 l GS for overall unrotated planes, the one loop Kahler potential involves 
all planes, whether rotated or not. 

2.2 Soft supersymmetry breaking terms 

The scalar potential of the N=l supergravity a— model is composed of F- and D- 
terms, V = Vp + Vd [40], defined as: 

V F = e K (D l W(K- 1 f J D 3 W - 3\W\ 2 ) = (F t G t] Fj - 3e G ), V D = ^ef/" 1 )^^, (7) 

where G = K + log \W\ 2 , D Z W = W % + K % W = d % W + Wd z K, D a = G^T^Aj, the 
moduli and matter fields are labeled by the index i and Fj = e G / 2 (G~ 1 ) i jGj identify with 
their auxiliary components. Following the hidden sector approach, we assume that the 
superpotential includes a non-perturbative part which lifts flat directions in S and Tj by 
inducing a non-trivial minimum of V characterized by non-vanishing VEVs for dilaton 
and moduli scalar S, Ti as well as auxiliary field components Fs,FT t - We shall follow 
the approach of Brignole, Ibanez and Muhoz [19] in which local supersymmetry breaking 
with < V >= is assumed to be saturated by the dilaton and moduli auxiliary fields. 
Generalizing the parametrization of ref.[19], we represent these quantities as follows: 



(G^) S§ F S = V3m g e ias sin#, (G^) T . T .F T . = V3m g e iaT eiCOS0, (8) 



under the restriction ^V e 2 = 1, necessary to ensure < V >= 0, and the simplifying as- 
sumption that the anisotropy parameters ei\i = 1,2,3] are real. The gravitino mass is 
m g = e ' 2 and 9 is interpreted as a dilaton-moduli mixing angle. The complex phases 
as, ax, along with possible non- vanishing imaginary parts of VEVs, ImS, ImTi, represent 
CP violation parameters. The isotropic case, which is considered in ref.[19], can be recov- 
ered by identifying the three moduli components, hence setting F^ = Ft, e* = 1/y/S. 

The case involving the one-loop corrected Kahler potential, eq.(5), can be treated in 
close analogy with the above tree level case by modifying the parametrization as follows 

[19]: 

i 
l( Fs+ *Y 6 JIs]^_) =v /3m e i°s sin ^ (i + ^h§Vjj±^ = ^ m jar e . 

Y\ {4x)2Z-'T i + T i ) g V (^) 2 YJ Ti + Ti g 

(9) 
The isotropic case results by setting, e* = -js,Sq S — 3^; $gs t+t- = ^ GS t+t • ^he 
parametrization (8) can also be generalized to allow for a matter component of the Gold- 
stino by introducing an additional mixing angle 9 A and writing [19]: 

(G^F)s =v / 3m g e MS sin#cos# A , (G* F) n = V3e iaT m g ei cos 9 cos 9 A , 
(G^F) A =V3m g e iaA sm9 A , 

the one loop case involving again a direct generalization of eq.(9). Auxiliary and scalar 
components are not independent, of course. They might be related by specifying the non- 
perturbative superpotential. Consider, for instance, the gaugino condensation approach, 
in which approximate elimination of the composite gauge fields associated to a subset 
l\ A , G' A of hidden gauge group factors, leads to the following toy model superpotential 
[13-15], 

W np (S t Tj = ^T/ue-^T J] [r,(Ti)] Ui , (11) 

A i 

where h A are constant parameters and v^ = —2 as required to obtain a modular invariant 
function, G np = K + \og \W np \ 2 , free of singularities in the Tj-planes. By virtue of modular 
invariance, the field space of the variables T can be restricted to the fundamental region, 
F = [\T\ > l,Im(T) < 1]. The relation between anisotropy parameters e$ and the ri 
parameters, defined by riTi(x) = T(x), can then be obtained by calculating the auxiliary 
fields associated with eq.(ll): 



{G*F) T . =m g [-l+u l (T l +T l ) — )J =m g 



00 r 

-l-^«(T i+ T 4 )fl-24^ U<h 



12 iy l '\ ^l~q\ 

(12) 
Since for Tj G F, Re(Ti) > \/3/2, and the moduli dependence involves the variable 

qi = e~ 27zT \ we see from the size of the exponential suppression, q^ < e _7rv ^ 3 rj 4.3310 -3 , 
that the relevant regime where the infinite sum in eq.(12) needs to be evaluated is that of 
large T$. Assuming that all Tj lie in the large compactification radius regime, one finds: 
— ~ — . For an isotropic orbifold, the restriction of equal r$ = r = 1, likewise implies equal 

ei = l/\/3. For an anisotropic orbifold, say, with T\ ^> T 2 ,T 3 , or r\ <C r 2 ,r 3 , assuming 
the large-T limit, leads to: ei ~^> e 2 , e 3 . 
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Let us now turn to the implications concerning the soft supersymmetry breaking 
interactions. We identify the string low energy observable sector with that of the minimal 
supersymmetric standard model and assume the presence of a perturbative superpotential 
describing the renormalizable trilinear and bilinear Yukawa couplings of quarks, leptons 
and Higgs bosons left chiral superfields, 



W = X^Q % eH 2 m + X d ij Q i {-e)H 1 D c j + X^L^-e)^ R c + pH x eH 2 



t-3 



ij 



(13) 



■ff°- 



such that Qi — Q'),Lj = (^),.ffi = ( 1 J),H 2 = (jo) denote quarks, leptons and Higgs 

weak doublets left-chirality fields, with the invariant combination, QeH = Qae^pHp, e a/ 3 
1] being the antisymmetric symbol; U c , D c , E c denote antiquarks and antilepton weak sin- 
glets left-chirality fields and the indices i,j = [1,2,3] label the squarks or sleptons gen- 
erations. We assume a one-to-one correspondence with massless string modes, so that 
modular symmetry implies that Yukawa coupling constants Af^T^), [x = -u, d, e] transform 

under SL(2, Z) m as modular forms of weight, (— 1 — n^ — n™- — n^ 2 ) for x = u and similar 

formulas for other cases. For a given orbifold, these are calculable functions of the moduli 
fields depending on characteristics of states, such as oscillator excitations, reference plane 
and, for twisted sector states, coordinates of fixeds points [39,41]. Adding the perturbative 
superpotential (13) to the non-perturbative one and substituting in the scalar potential, 
eq.(7), yields a tree level soft supersymmetry breaking potential of the familiar form, 



ei2 



V soft 



z,2 



If) 



A a Ai + -M a X a X a + 



m. 



A^Q t €H 2 m + A'XgQii-e^D 



+A e lj X l e J L i (-e)H 1 E^ + B^H x eH 2 + c.c. 



(14) 



with the following formulas for the soft breaking coupling contants in the tree level case [19], 



M, 



rrir 



° Re(fa) 



V3e~ ias sme(^)Re(S) - T^e~ iaT cos6 ^KM^C^T,] 



dS 



m 2 =m 2 Jl + n a " cos 2 6>), 



A 



ctpj 



B, 



V3e~ tas sin 9 - e~ iaT cos 9 1 e' + n' a + n'p + ri^ - v^ ^(T* + T^e* 



xv 



ct{3~f 



IV 



aji^j 



1 - V3e~ tas sin 9 1 - ^-(S + S) 



Vs 



//, 



- e~ iaT cos 



e(e' +ri Hl +n' H2 - y/3^2 ei(Ti +T,) — )• 



(15) 
The index a attached to low energy matter fields, denoted generically as A a (x), subsumes 
gauge and generation quantum numbers and the index a = 3, 2, 1, attached to gaugino 
fields, X a (x), labels the standard model gauge group factors SU(3) x 577(2) x £7(1), in 

the indicated order. We have introduced the following notations: n Q = 3 J2i e l n a^ n a = 
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\/3E^?C C = v / 3E i e / i ,G 2 (T) = G' 2 (r)-^, G 2 (T) = -4tt^ « tt 2 (1/3 - 
8 E n ^g n /(l— Q n )) ; ws = dw/dS, with u; = w a p 1 A a ApA 1 a generic notation for the super- 
potential terms in eq.(13). The isotropic orbifold case, which is treated in ref.[19], is recov- 
ered by setting, a = 1/V3, which gives: e' = 3, n a = n' a = n a = Ej n ^ v^ E; e i( T i + 
Ti]JTi = (T + T)fr. The one loop improved case is described by identical formulas to 
eqs.(15) with the substitution: 



Cv 



e t 



i 

1 _l_ GS 

1 T (47r) 2 r 



The isotropic one loop improved case is recovered by performing the substitutions described 
above for the tree level case, together with the formal replacement [19]: 

Oc* c l 

cos6> — > cos6*(l H )~2, sin 6* — > sin 9. 

2Ati 2 Y 

2.3 Renormalization group analysis 

Distinct physical theories are parametrized by renormalization group trajectories 
which trace the scale dependence of the gauge, Yukawa and soft supersymmetry breaking 
coupling constants . These are described by first-order differential equations in the scale 

variable, t = log -y?f , with Mx the string unification mass scale and Q a floating renor- 
malization scale, so that a solution to these equations is fixed uniquely once boundary 
conditions are specified at Mx, as for the quantities quoted in eqs.(15), or at some other 
scale. While results for the renormalization group equations are available in the litera- 
ture up to two-loop order, in the present work, we shall restrict ourselves to the one-loop 
equations, and rely on refs. [21,23,42] as our main sources regarding notational conventions. 
The one loop gauge coupling constants and gaugino masses are then decoupled from the 
other parameters. Their scale dependence can be expressed explicitly, 

n, m - 9l(t) q a (0) M a (0) 

(Ait) 2 1- f3 a t + a a {0)A a 1 - (j a t 



where: 



P a = b a a a (O) = b -^0^, A a = 6 a + '£b'*]Qg[{T i + T i )\T,{T i )\*]. (17) 



The string threshold corrections A a include a moduli-independent part 5 a , to be neglected 
in the following, along with the familiar moduli-dependent part [9]. The /3-function slope 

parameters b a are defined as: dg a /dlogQ = —b a g^/(4:TT) 2 , the values for the standard 
model group factors being: b a = (3,-1,-11), subject to the tree level normalization 
convention, gl(0)k a = g\- 

The renormalization group equations for the Yukawa coupling constants involve the 
gauge coupling constants and a dependence on generations through the Yukawa couplings. 
Motivated by the large hierarchy between third and first or second generations, we shall 
adopt here the approximation where all entries in Af-, [x = u, d, e] are neglected relative 

to the (3, 3) entry. Denoting the corresponding (33) Yukawa matrix elements as A t , A&, A T , 
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the index standing for top- and bottom-quark and tau-lepton, respectively, one obtains 
the simplified third generation evolution equations [23]: 

^^ =Y t (t) (y ua a a {t) - 6Y t (t) - Y b (t) 

dY b (t) 
dt 

dY T (t) 
dt 

where Y x {t) = $$,[x = t,b,r] and y ua = (f , 3, f ), y da = (f,3,|),y e a = (0,3,3). An 
analytic solution to these equations exists in the approximation X t ^> A;,, A T [23], 

Y t (t)= m^l , Y b (t)= fb{ ° )m) 1; Y T ( t) =Y T (0)E 3 ( t) , 
tK) l + 6Y t (0)F 1 (ty bK) (l + 6Y t (0)F 1 (t))r — 



--Y b (t) (y da a a {t) - 6Y b (t) - Y t (t) - Y T (t)\ , (18) 

--Y T (t) (y ea a a (t) - AY T {t) - SY b (t) 
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where: 



Via , 



E 2 (i) = Ei(i)(l-ft<)™ 7 , E 3 (*) = (l -&*)"(! -ft*)" (W) 

The supersymmetry breaking parameters mix together as well as with the gauge and 
Yukawa parameters. For the Afj,[x = u, d, e], we shall adopt an analogous approximation 

to the above, namely, dropping all but the (i,j) = (3,3) entries. This approximation 
is justified by the fact that these quantities always enter through products, AfjXfj. The 

simplified evolution equations for the (3, 3) parameters, designated as A t , A b , A T , read [23]: 

d -^- =y U aa a (t)^ _ 6Y t (t)A t (t) - Y b (t)A b (t), 
at m g 

^^ =y da a a (t)^^- - 6Y b (t)A b (t) - Y t (t)A t (t) - Y r (t)A T (t), (20) 

at m g 

^7^ =y ea aa(t)^l _ 3Y b (t)A b (t) - 4Y T (t)A T (t). 
at m g 

None of the other matrix elements, Afj, [x = u,d,e] with (i,j) ^ (3,3), are coupled to 

the Af 3 (t), so that these parameters have no other renormalization scale dependence in 
the present approximation except for that arising from gaugino loops. This affects the 
diagonal elements by a generation independent contribution analogous to that arising in 
eqs.(20). The approximate solution can be obtained by a simple quadrature, 

A*.(t) = A* j (0) + yxa a a (0) Ma{ ° ) * Stj [x = u,d,e; (i,j) # (3,3)]. 

trig J- Hat 
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The renormalization group equations for the scalars masses include contributions 
from gauginos loops as well those from scalars and Higgs bosons loops. In the above 
approximation of dominant A| 3 (t) Yukawa couplings, the scalars mass matrix elements, 
(i,j) = (3, 3), and (i, 3), (3, i), [i = 1, 2] decouple from other Yukawa matrix elements and 
have a scale evolution governed for the (3,3) entries by the equations [21]: 

cl \ m Q)33 _ 2 t> / ~ 2 , ~2 i -2 , 2 a2\ t> / ~ 2 , ~2 , -2 , 2/|2\ 

^ =z Qa a a M a -y t (m Q + m Q c +^ 2 +m ff AJ - Y 6 (m Q + ra 6 c + Mi +m g A b ), 

mQc --z Qca a a M 2 a - 2Y Q (m 2 Q + m 2 Qc + (^\ + m 2 A 2 Q ), 



dt 



— k =z La a a M 2 - Y T {m\ + m 2 rC + Jl\ + m 2 g A 2 T ), 
-^p =z LCa a a M 2 - 2Y T {m 2 T + m 2 rC + Ji\ + m 2 g A 2 T ) 

and for the entries (i,j) with [i = 1, 2, j = 3] or [i = 3, j — 1, 2] by the equations, 

d 



-fA m Q,Q c ,L,L- 



dt 



iij 



^(Y t + Y b ),Y Q ,^Y T ,Y T 



(m c 



,Q c ,L,L c )ij- 



(21) 



(22) 



We have omitted writing the argument t for the various coupling constants and have 
denoted Q = [t, 6], Q c = [t c , b c ], L = [r, v T ], L c = t c with zq u = 4C2(-R Q ), corresponding to 
the quadratic Casimir operator in representation R a for the standard model gauge group 
factor a = 3,2,1, such that z Qa = (^,3, ±), z u c a = (^, 0, ^), z d . a = (™,0,%),z L c a = 
(0,3,l),z eCa = (0,0,4). 

For the scalars masses matrix elements, (i,j) with [i = 1, 2, j = 1, 2], we shall include 
the contributions of the Yukawa interactions perturbatively by retaining the terms with 

leading power of t = log-g^. The gauge interactions contributions, being generation 

independent, have the same form as those appearing in eq.(21), and are then included 
exactly. The relevant formulas read [21]: 



(rfiQ(t)h « (m 2 Q (0)) lJ + 5 lJ G Qa (t) + t 



[~ ^({K^rh 2 Q }).. - (Km 2 Q \l 



uJtJ 



-l£(\ u ^l)ij - m2 g[ (X U A U )(X U A U ) J( I ] + [X u ->• A d , [i2 ->• A*i] 



IJ 



{m Qc {t))ij « (m Q (0)) u - + 8ijG Q c a (t) +t 



({X^XQ.mlo}).. 



-2(X^ Q m 2 Qc X Q ) tJ -2[^y^ Q X Q ) ij - 2m 2 g [{X Q A Q )\{X Q A Q ) 



>3 



(23) 



where we understand that terms on the right hand sides are evaluated at the electroweak 
scale and where 



1 



G Qa it) = -z Qa a a (0)M^(0)f a (t), [f a (t) 



(2-P a t)t. 

(i-M 2 - 
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stands for gaugino loop contributions. The slepton mass matrices fh 2 L ,m 2 Lc are obtained 
by analogous formulas to those for mQ,ifiQ C above by substituting slepton mass matrix 
for squark mass matrices and letting A^ — » A e , X u — ► 0. 

The mass parameters ft 2 = fi 2 — [i 2 , appearing in formulas above for scalar masses, 
refer to the coupling constants in the Higgs bosons scalar potential, 

V(H t ) = ^l 2 + £\H 2 \ 2 - fiiHtH* + c.c.) + ^^(|^| 2 - |# 2 | 2 ) 2 , (24) 

o 

where fx 2 = m 2 H . + fx 2 [i = 1,2], fx 2 = —B^m g ^. For completeness, we quote the scale 
evolution equations for these Higgs sector parameters, again in the simplified version in 
which our calculations are ferformed, namely, predominant (3, 3) trilinear couplings [23,42]: 



dy 2 (t) _ 2 
dt _/X 
dji\{t) 



3d 2 + di - 3{Y t + Y b + -Y T 



dt 

di4(i) 

dt 

d/4(t) 



--3a 2 M% + a x M\ + ( 3d 2 + &i - 3(Y t + Y b + ^Y r ) ) /j 2 - 3Y" b (m 2 + m\ c + [i\ + m 2 g A 2 b ), 



=3« 2 M| + a x M\ + ( 3d 2 + &! - 3(Y t + Y b + ^Y T ) U 2 - 3Y t (m 2 + m 2 c + Jx\ + m 2 g A 2 ), 

1 



I- 
3 

I- 
3 



- (3d 2 M 2 + a 1 M l - 3m g A t Y t )ii + ^(3d 2 + &i - 3Y t - 3Y b - Y T )^ 2 . 

(25) 
Let us recall here that the weak hypercharge D-terms add an extra contribution in the 
scale evolution of scalars masses [22] , 



nun 

^ = -Y a a 1 (t)S(t) + - 
dt 



S(t) 



g(g) 

I -Pit 



S(t) = 2J(toq - 2m 2 JC + m 2 D c - m 2 L + m 2 Ec ) - m 2 Hl + m 2 H2 = Trace(m 2 a Y a ), 

gen 

where Y a = (^, — 3, ^, — ^, 1, — ^; 2) * s ^ ne hypercharge of the low energy fields in the 

order 

(Q, U c , D c , L, E c , Hi, H 2 ), and dots refer to gauge and Yukawa terms quoted above. This 

term contributes only in cases with non- universal boundary conditions [22]. Although it 

may significantly influence the spectra, it was found in ref.[33] to affect negligibly flavor 

changing observables and so will be neglected in the sequel. 

At the electroweak breaking scale, which we shall identify with the Z-boson mass, 

Q ~ mz, tz = 2 log — — , scalars masses obtain additional contributions from Yukawa 

couplings F-terms and gauge couplings D-terms. Let us adopt the familiar notation in 
which the three generations of left chiral squarks (L) and anti-squarks (R) are combined 
together into a single six-dimensional coloumn vector, so that the scalars mass terms in 
the Lagrangian appear in the 6x6 matrix form, 



E (Ql>Q 



Q=U,D 



M 



Q2 

LL 



M 



Q1 

LR 



AfQzt m® 2 

1V1 LR 1V1 RR 



Q 



R 



+ IQ-L], 



(26) 
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The full squarks and sleptons masses, comprising the above supersymmetry breaking 
contributions evolved from the the string scale down to the electroweak scale, along with 
those induced at the electroweak scale, read [42]: 

(M?lh(t z ) ={rh%Mt z ) + {M Q M ] Q ) l3 T m 2 z cos 2/3 (^ - ^Q sin 2 w )8 ij , 
(Mf R ) l3 {t z ) ={m 2 Qc ) i3 {t z ) + (M^ Q M Q ) lJ Tm 2 z cos2(3^(^j sin 2 6 w 5 tJ ; 



1 , _/l-2sin 2 0i 



(Mtlhitz) =(fhlh(tz) + (M L Ml) tj ± -m|cos2/3 \ * W )6 i3 , 



2 * ' V I 



(M^ijitz) =(mic)»j(tz) + (MlM L ) l0 - m z cos 2/3 sin 2 W <5. 



y 



(^S)ti = [m^ + ^](M u )y, (M& e)2 )y = [m s ^' e) +^tan/3] (M (d>e) ) y ; (27.6) 

the upper and lower positions in eqs.(27.a) referring to the up-squark and down-squark 
cases, respectively, or to slepton and sneutrino cases, respectively. We denote by M x , [x = 
u, d, e] the quark and lepton mass matrices, M u = \ u (tz)^%, ^d,e — ^d,ei^z)^hi where 

< Hf >= Vi/\/2, such that v 2 + v\ = v 2 , tan/3 = v 2 /vi, m 2 z = (g 2 + g 2 )v 2 /4. 

The condition of radiative electroweak symmetry breaking, as derived by minimizing 
the Higgs bosons potential, eq.(24), involves solving the two equations [23,42]: 

l 4(t z )+m%l2' A{tz)+A(tz) 

For completeness, we also quote the mass terms of charginos and neutralinos, L mass = 
-[ X - ] M cX + + c.c] - ±x 0T M nX °, with 



■i 



givi 


91 V2 v 

2 \ 


2 


52^1 


32 v 2 


2 


2 





fi 




J 



/Mi 

/ M 2 ^a \ f M , 

M = I ^ I M = 

jw c I g 2 vi _,, I 5 lvl n 

V V2 ^ / I 

referring to the basis, x ± — (—iW ± ,H2i), X° — ( — iB°, — iW°, H®, H®)- The contribu- 
tions of super-box or super-penguin diagrams to flavor changing neutral current processes, 
with external quarks and leptons, are most conveniently described (in the so-called inser- 
tion approximation) in terms of generation mixing terms for scalars mass matrices in the 
following representation: 

M 2 N = V M M 2 MN V^ [M, N = L,R] 

corresponding to the super-CKM mass bases for the quarks and leptons superfields with 
generation diagonal .D-terms Yukawa interactions, Q*T a \ a qj5ij + L*T a \ a ej5ij +c.c The 
transformation matrices V^ [M = L,R;q = u, d, e], relating current and mass bases for 
quarks and leptons, are defined in the conventional way so as to satisfy, M q = Vj* {M q ) ( nagV R [q 
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u, d, e] . The above notation emphasizes the important fact that flavor changing observ- 
ables are sensitive to supersymmetry breaking as well as to Yukawa interactions. Note 
that only the relative product, identified with the (Cabibbo-Kobayashi-Maskawa) CKM 

matrix V = Vqkm — V^V L , is physically observable if one restricts to fermions weak 
interactions, while additional information on the structure of the fermions mass matrices 
is needed to describe the physical flavor changing paramaters. 

2.4 Yukawa coupling matrices 

To obtain boundary conditions for the Yukawa coupling constants at the string unifi- 
cation scale, one might use, as in the so-called bottom-up approach, the available exper- 
imental information for quarks and leptons masses and weak mixing angles to determine 
the Yukawa coupling matrices at the electroweak scale and next evolve these up to the 
string unification scale by integrating the renormalization group equations [21]. Here, we 
shall follow a simpler approximate procedure: We assume some fixed choice for the Yukawa 
coupling matrices at the unification scale, with a structure consistent with orbifolds pre- 
dictions [39,43], and deduce from this the implied form of V^ at electroweak scale. The 
analysis of Casas and Munoz [43] of the generation dependence of Yukawa matrices for 
orbifolds favors the following structure, 



A 9 (M 



x) 






a 





a 


A 


C q 





Cn 


B 



[q = u, d, e] 



(29) 



the entries here being real numbers, a choice which can always be achieved by performing 
suitable chiral transformation field redefintions. The choice of eq.(29) has a suppressed 
value for the first generation, assumed to arise from non-renormalizable interactions, and 
a nearly diagonal structure for the second and third generations couplings, assumed to 
arise from renormalizable couplings, the suppression of off-diagonal terms being caused by 
a proper choice of widely spaced fixed points. The assumed hierarchy of parameters in 
eq.(29) is: B ^> A ^> a, c q . For top-quark masses, m t < 100 GeV, the scale evolution of the 
Yukawa matrix is dominated by gauge bosons loops and can be represented approximately 
as an overall scaling factor determined by the gauge coupling constants, 



M q (t) = F q (t)M q (0), F q {t) 



nl &a\t) \y a /b a 



Full formulas for the scale evolution factor are quoted in ref.[43]. For larger top-quark 
masses, the top-quark Yukawa couplings compete with gauge interactions and affect mainly 
the matrix elements, Af- with i — 3 or j — 3, corresponding here to the parameters B(t) 

and c q (t). Considering, for definiteness, the down-quarks case, the mass eigenvalues for 

a Yukawa matrix of the form of eq.(29) are: rrid ~ a 2 /A,m s ~ A, nib ~ B, so that 
only nib is scale-dependent in the approximation of dropping Yukawa interactions in the 
scale evolution, while c q is scale dependent if one includes third generation couplings. For 
symmetric mass matrices, as is the case here, left and right transformation matrices are 
equal, V£ = V£, and given by: 



v£ 



V R 



( 



1 



I raq | 



^/m d m s c d 
5 

rut 

b 



Crf 

m b 



^m d m a c d \ 
m b m s 



Cd 

m b 



(30) 
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so that the transformation matrices are fully specified, up to a single unknown param- 
eter, Cd = Cd(tz) [43]. Note that (Mdjdiag = (-mrf,m S) m6), and that the negative 
mass eigenvalue can always be absorbed, if so desired, by a chiral transformation, say, 
d,L — > —d,Ljdji — > d,R. We have only quoted above explicit formulas for the down-quarks 
case. We shall assume that similar structures hold for up-quarks and leptons mass matri- 
ces, so one obtains for these cases analogous formulas by replacing (d, s, b) — > (u, c, t) or 
(e, /j,, t) with two additional free parameters, c u (t), c e (t). The values of these parameters at 
the electroweak scale can be obtained by identifying V^V L with the CKM-matrix, which 
supplies four relations expressing c u (t),Cd(t) as linear combinations of the CKM matrix 
elements, V u b, V c b or Vta, Vts- We shall use in our calculations the above form of Yukawa 
couplings for up and down quarks and for leptons, and following ref.[43], tentatively set the 

c q (t) parameters as follows: c u (tz) ~ rn s \p2, Cd(tz) ~ fn Cl ci(tz) = m^. Of course, other 
choices for the Yukawa matrices will lead to alternative expressions for the transformation 
matrices. For symmetric Yukawa matrices, under restrictive assumptions on the number 
of texture zeros, the number of choices is known to be limited to only five cases [44] . 

3. Results and discussion 

In this section we shall present numerical results for the dimensionless parameters 
defined by: 

iCf'2 yQ f~rq2 yq\ 
X 1V1 MN V M 1V1 MN V N /oi\ 

<JMN — — ^T~ — z-s , {oi.) 

where the mass factor in the denominator in represents some average scalar superpartner 
mass, which we shall set for the equal chirality cases 8mm to the weighted trace of the 
mass submatrix of fixed chirality M = L, R and for the mixed chirality case 5lr to the 
weighted average of the full mass matrix trace. 

3.1 Inputs 

The basic parameters entering the calculations are set as follows: 

(a) Masses: m z = 91.17GeV,v = 246GeV;m t = 170GeV,m b = 5.6GeV,m T = 
1.784GeV, m c = 1.35GeV, m s = 0.175GeV, m M = 0.1056GeV, m u = b.lMeV, m d = 8.9MeV, m e 
0.51MeV. 

(b) Gauge coupling constants : g\{mz) = 0.127, g\{mz) = 0.425, g\{vaz) = 1.44, gx = 
A=,M X = 5.27g x l0 17 GeV = 3.7310 17 Gey, with t z = 2 log ^ = 71.9. The boundary 
conditions a a (0) are evaluated via eq.(16). 

The dilaton VEV is set at < S >= -\- = 2. The numerical values for the moduli 

VEVs will be chosen in the various cases to be considered below on the basis of fits to 
gauge coupling constants unification, and specifically to a s (tz) and sin 2 0\y(tz) [8,45]. 
Recall that the radiative corrections to these quantities involve the linear combinations, 
A 'i = fci^i _ b '^ B 'i = b 'i + h 'i _ kij^b'j, so defined as to be free of the 5 l GS . Note 

that the typical values of order T ps 10 that are found in these fits [8] lie an order of 
magnitude above those found in mimimizing the scalar potential in gaugino condensation 
models [17]. In the present work we shall not consider CP violation effects and therefore 
will set to zero the imaginary parts of dilaton and moduli VEVs and phases of auxiliary 
fields, as = cut = 0. We also shall neglect contributions to supersymmetry breaking 
parameters, eq.(15), involving derivatives with respect to dilaton and moduli fields of the 
Yukawa coupling constants. 

The Green-Schwarz parameters S l GS , representing gauge group independent contribu- 
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tions to duality anomalies cancellation, are calculable for orbifolds on a model-by-model 
basis via the /3-function slopes of N = 2 suborbifolds with z-plane fixed. The typical values 
in a minimal standard model orbifold can be calculated from formulas given at the end of 
Subsection 2.1. One finds for the diagonal sums, Sqs = X^ ^GS' values °f positive sign in- 
side the interval, 5qs ~ [5 — 10]. (Note that our conventions for the parameters b a , b*, 5 l GS 
differ by a sign from those of ref.[19].) This contrasts with the range 5qs ~ [25 — 50], 
favored by gaugino condensation models [17]. 

The parameters h* also are model dependent. If large moduli VEVs are used, then, 

because of the exponential dependence on Tj, this strongly constrains the b* of overall 
unrotated planes if one is to avoid excessive threshold corrections. Indeed, expressing 
these corrections to gauge coupling constants in terms of an effective string unification 



scale M' x , one finds: 






M' x 


-u 


(T, + T0l^m)| 4 



Y[(Ti + T^ie < , (32) 



where the large radius limit was used in the second step. We see that a reduction of 

order 10 in the string scale requires for the exponent — ^f- a positive value of absolute 
magnitude close to unity. An analogous constraint applies for the more relevant exponent, 



~t - 



~ 2(6 k -b k) associated to the unification of two group factors, G a and G\,. 

The experimental information for 5mn can be obtained from flavor changing neutral 
current observables involving mass differences in the neutral mesons systems (K — K,B — 
B, D — D) and branching ratios of transitions such as qi — > qj + 7 or 1$ — > lj + 7, [i^ j]. 
Upper bounds are derived by assuming that these observables are solely accounted in 
terms of contributions of superbox or superpenguin diagrams with gluino and squark, or 
photino and slepton, exchanges. Denoting the average superpartner scalar mass by m, the 
gluino and minimal neutralino masses by M3 or My, then the calculated mass differences 
5m k, 8m b, Sm^ and branching ratios B have a dependence on superpartners masses which 

M 2 

can be separated into two variables: m and xs n = -^ L . A study of formulas quoted in 

refs. [20,34] reveals that 5rriK,B oc -^t-Fb(x), with Fb(x) smooth functions of x, so that 
the box diagrams bounds imply the scaling law, 5m n oc m, [M, N = L,R]. Using x = 1 
as a reference value, then account of the x— dependence results in stricter (weaker) bounds 
depending oni < l(x > 1). The penguin diagrams contributions to branching ratios, 
neglecting the contributions of A— and fx— coupling terms [20], read: B oc 5 MM Fpi(x)/m 4 
or 8 2 LR Fp2(x) 1 'm 2 ', which implies the scaling laws, 5mm oc m 2 , 5lr oc m, times smooth 
functions of x. 

The numerical values of presently known bounds for 5m n , as obtained from results in 
refs. [20,27,34], are quoted in Table 1. More stringent bounds are found for the quantities 

in Table 1 designated as, 5ij = (5ll5rr)^ , which correspond to contributions implying 
simultaneously non-vanishing chiral left and right mass terms, 5ll5rr 7^ 0. A similar 
situation holds for the case of simultaneously non-vanishing (Slr)^ and {5lr)ji matrices. 
The bounds in Table 1 for 5lr correspond precisely to the geometric average of these 
matrix elements. 

3.2 General features of model 
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It is useful to identify the various sources of off-diagonal contributions to the matri- 
ces (Smm)ij- The F- and D-terms (second and third terms in the formulas of eq.(27.a)) 
yield purely diagonal contributions. The solutions of the renormalization group equations 
(first terms in eq.(27.a)) give contributions the scalars mass matrices with the approximate 

structure, (M^)y(*) « (m Q %(0) + G Q (t z )5 tJ + T, q ^u,dvf{ M i M l'^^ under the sim- 
plifying assumption that the contributions from boundary conditions, gauginos couplings 
and Yukawa couplings add up linearly. Off-diagonal matrix elements cannot arise from the 
family independent gauginos couplings. They also cannot arise from the boundary con- 
ditions mass matrices if these feature family universality (to 2 - is multiple of identity) or 

alignment with the fermions mass matrices (VM'fh 2 vIN (0)V^ [ are diagonal matrices). As is 
also well known [20] , the radiatively induced off-diagonal terms from Yukawa couplings are 
numerically significant only for the sub matrix Mf^ L , yielding for flavor changing parameters 
the approximate result: 



( d ) / \ 2 

(#L)y «|rkM*a,v » y^v^Vsj, (33) 

" l d V / a rn d 



m d V 

id) 



with yu of order unity and V = Vckm- (The chirality L is singled out because only 
left-chirality scalars of given charge interact with scalars of both charges and the charge 
d because of the predominant top-quark Yukawa coupling.) For these flavor changing 
terms of standard model origin one expects then important suppression due to small 

CKM mixing angles. Using numerical estimates quoted in ref.[34] for yu , gives: 5f 2 ~ 



7.3 x 10~ 4 , Sf 3 «1.4x 10" 3 . 

For the mixed chirality parameters, the relevant terms in the interaction basis which 
contribute to off-diagonal elements of (olr)ij can be written approximately as: 

■m g cos9(m+nj)+M a (0)y xa a a (0)Sij 



1-Pat 



"(34) 

This formula clearly illustrates the fact that non-universal contributions (first term) could 
be masked by gaugino contributions (second term), in the presence of large gauginos masses. 
Because of the dominant (3, 3) entries in the fermions mass matrices, the expected pattern 
is: (SIr)23 > (<$£fl)l3 > ($lr)i2- 

Returning now to the unmixed chirality case, a simple estimate of the generation 
non-universality of boundary conditions can be performed by including the contributions 
of gaugino loops in the scale evolution only, while dropping entirely those of Yukawa 
couplings. The physical mass matrices in the weak interaction basis are then diagonal 
ones and given by: 

(M x2 )ij (tz) =m 2 Jl + n xi cos 2 9 + ^ sin 2 9 ^ z xa a a (0)f a (t z )) <%, [x = Q, e] (35) 

using self-explanatory notations, with f a {t) defined at eq.(23) and dropping momentarily 
the suffix MM = LL, RR. Performing the transformation to the super-CKM basis under 

the simplifying assumption of negligible off-diagonal matrix elements M 2 - ps [i ^ j], 

one obtains the approximate formulas: 

M t f = {V^M x2 vS)ij « Mr 2 V^ J +MffV^ J HM^ 2 -M^)V^ 3 V^ 3j +- • • , [x = Q,e] 
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where dots refer to contributions proportional to (Mf? — Mff)- Given the structure of 
Vm, eq.(30), the expected pattern for the flavor changing contributions associated with 
the boundary conditions is: 5\2 ~ ^23 > ^13 • Neglecting the F- and D-terms, one obtains 
for (i,j) = (1, 2) or (1, 3) the rough estimate: 

s* ~ 2 (v x h M ^ - M % ~ cos2 - (v^) Uxi - Uxj (36) 

ij [ M)lJ M** + Mff 19 ^ M N+^(^±^-18)cos 2 0' W 

which shows that gauginos can have a strong dilution effect on non-universality at the 
unification scale, to the extent that the mixing angles 9 lie not too close to or n. 

We shall now focus on the right chirality d-squarks case since this possesses simplifying 
features allowing to find upper bounds on the flavour changing parameters. Only the diag- 
onal matrix elements here get renormalized by gauginos loops. Since this is proportional 
to the identity matrix it follows that the right chirality d-squark mass matrix has non- 
renormalized non-diagonal terms. Let us write down this matrix at the compactification 
scale in the basis where the d-quark Yukawa couplings are diagonal: 

MS&j = ^ (S tj + cotan 2 # J2 VLfa + 1)V&) , (37) 

where mo = M a (0) is the common gaugino mass at unification scale and we denote, for 
notational convenience, rii = udi- Defining the renormalized average mass by means of a 
numerical estimate for the gaugino loops contributions, 



,2 



m 0/io , n l + n 2 + n 3 + 3 



m av = ^"( 18 H § COtan e )' ( 38 ) 



we find the flavor changing parameters, 

(fid ) _ Efc Vik n kV k j ,, q v 

This ratio explicitly depends on the arrangement of modular weights. Using the ansatz for 
the V d matrix, eq.(30), the numerator can be evaluated. As we are looking for an upper 
bound we choose the arrangement generating the largest contribution. We also safely 
neglect the third family angle compared to the Cabibbo angle. Using unitarity the result 
can solely be expressed as a function of modular weight differences. Moreover the 9 angle 
can be constrained by requiring that the scalar masses are positive at the compactification 
level. This implies that tan 2 6 > — fT7 _ 1 _ l _ n , where n m is the lowest modular weight. This 
yields the upper bounds: 

(Srr)i2 < 



(n m + l) ' 



m d Arai 



m s -19(n m + l) + Ani + An2 



d V2m d m s An2 

(<Wl3 - m b _i 9(nm + 1) + An^ ( 4 °) 

f X d \ / ft™* An 2 

{0 RR )23 < V2 



m b -19(n m + 1) + Ani + An 
20 



_> 



where Am is the modular weight difference involving m, e.g., Ani = n\ — 713 if n m = 773, 
and similarly for A772. Notice that the order of magnitude of these mass insertions is given 
by the ratio of d-quark masses. Effectively these ratios are not very different from the ele- 
ments of the CKM matrix. Indeed the ratios of the u-quark masses composing the u-quark 
Yukawa coupling diagonalising matrix are negligible compared with their d-quark counter- 
parts. Thence (S RR )i 2 is proportional to a very good approximation to sin 6>c- Numerical 
values can be extracted from eq.(40). For instance, using m = —2, 772 = —3, 773 = — 1, 
one obtains (S RR ) 12 < 6.10" 3 , (S RR ) 13 < 2.5 10" 4 , (S RR ) 23 < 1.1 10" 3 . These analytic 
bounds are modified numerically by the subleading Yukawa couplings. Nevertheless their 
orders of magnitude is retrieved. 

3.3 Non-universal modular weights assignments 

Most examples of solutions quoted by Ibanez and Lust [8] are generation indepen- 
dent. It is clear, however, that relaxing the generation independence constraint could only 
increase the number of solutions. Rather than repeating the type of analysis in ref. [8] to 
seek for generic features of solutions involving a generation dependence, we choose here to 
consider three representative solutions, taken among the examples quoted in ref. [8], and 
consider generation dependent discrete perturbations of modular weights by a few units. 
We do not require consistency with anomalies cancellation and unification constraints. In 
all cases discussed below, we set the standard model Kac- Moody levels at: k a = (1, 1, |). 
The following three cases will be considered: 

(1) Case I: (ref. [45] and ref. [8], eq.(4.6)) This refers to a general, unspecified orbifold 
case with one modulus, constrained by gauge unification only and by a restricted interval 
for modular weights : —3 < n a < — 1. It can be realized in terms of an isotropic orbifold or 
an anisotropic one with T\ ^>T 2 ,T 3 , both cases with all planes overall unrotated, N unr = 
3, or with two rotated planes. The generation independent isotropic case, disregarding 
therefore the anomalies cancellation condition, admits an unique solution consistent with 
the unification conditions on the strong and weak angle coupling constants [45]: uq = 
no = —l,n[/ = — 2,77l = he = —3,77^ + uh 2 — — 5, — 4, with Re(T) ~ 16. (The 
suffices U,D,E here obviously refer to the corresponding antiparticles.) We choose here, 
tih 1 = — c 2,TiH 2 = —3. Using the above values of modular weights yields: b' a = (6,8,18). 
The exponential dependence of threshold corrections and large VEV of the overall modulus 

entail a careful tuning of Sqs- We shall set 5qs — 7, i.e., 6 a * = (— 1, 1, 4^, which gives 
j-^- ~ (10 _1 , 10~ 3 , 10~ 2 ) for [a = 3,2, 1]. We shall implement generation dependence in 

terms of two subcases: (i) Case I-a: n a i = n a — l,riLi — ul + l,^si = tle + 1, [i — 1, 2] 
and 77q,3 = 77 a , ?7l3 = 77_l, 77^3 = Tie with a = Q, U c : , D c ; (ii) Case I-b: Same assignments 
with Re(T) = 1.3. 

(2) Case II: (ref. [8], Section 4) Z' 8 orbifold with two rotated planes, % — 2,3, consis- 
tent with unification, Re(Ti) rs 24 and anomalies cancellation, 6 a 2 ' 3 = 0. This admits 
several generation dependent solution, of which an example is: 77q 12 3 = (0, — 1, 0), n l D123 = 
(-l,0,0),77[ ri = -i(0,l,l),n^ 2 3 = -|(6,15,3),^ 1 = -|(14,3,7),^ 23 = -i(14,7,3),n^ 1 
-(l,0,0),77 l E23 = -i(6,15,3),77^ = -i(4,3,3),77^ 2 = -|(14,3,7). The /^-functions 
slopes are calculated to be: b* — (| , | , ^), b^ — k a 5 l GS = |(_ 7 ), for i = ( 3 ). We shall 
choose 5q S = |, which gives Sqs — 8. The anisotropic torus with Ti ^> T 2 ,T 3 , e 2 ,es <^i 

ei ~ 1 has the modified weights, n'^ ~ V3n' a , defined at eq.(15), given as: 77g 123 = 

0' n D123 = — 3, ?%i23 = (0, —4, —4), n L12Z = ~~' n El 2 3 = ( — ^, ~j, —4),"-^ = ~ 2 , ^h 2 = 

— -£. We see that these effective weights cover a wider range than in the isotropic case. 
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We consider two versions based on this example: (i) Case Il-a: Isotropic torus with overall 
weights: n Q12 3 = n D12 3 = -l,nui = -l,nu2,3 = -3, n L i23 = -3,n E i = -l,n E 2,3 = 
—3,n//! = — 2,njy 2 = —3. (ii) Case Il-b: Anisotropic torus with T\ ^> T2,Ts, e2,e3 <C 
e\ ~ 1. To simplify calculations we represent this case by using the usual weights, but with 
values covering a wide interval of the type described above. We use: uq\23 = noi23 = 
n L i23 = (0, -2, -3), nc/123 = n E12 3 = (-1, -3, -S),n Hl = -2, n H2 = -3. 

(3) Case III: (ref.[8], Section 4) Z' 6 orbifold with one rotated plane, i — 3, consistent 
with unification, Re(T\) = Re{T2) ~ 10, and anomalies cancellation. The first condition 
yields ^2 i=1 2 & a * = (5>4, — ^-) and the second yields 6 a 3 = k a SQ S = —2. To achieve 
satisfactory gauge coupling constants unification we shall choose: ^2 i=1 2 ^gs = 6. Starting 
with the universal solution of ref.[8], uq = ud = njj = ul = n E = — 1, uh x = nn 2 — — 1> 
generation non-universality is implemented by considering two subcases, (i) Case Ill-a: 
Re(Ti) = 10, n a i = n a ,n a 2 = n a — l,n a 3 = n a — 2; (ii) Case Ill-b: n a \2 = — l,n a 3 = 
-2, a=[Q,U c ,D c ,L,E% 

3.4 Electroweak breaking, soft parameters and low energy spectrum 

We choose the dilaton-moduli mixing angle 9 and gravitino mass m g as our sole free pa- 
rameters and determine the remaining parameters, namely, tan/3 and //, by solving eqs.(28) 
at the electroweak scale with the solutions of the renormalization group equations for 
//|(t), [i = 1,2,3] obtained by integrating the coupled equations (20), (21), (25), using 
eq.(19) for Yukawa couplings. We shall present detailed results for the following choice 
of supersymmetry breaking free parameters: m g = 200 GeV, 9 G [0, n] and for Yukawa 
coupling constant : Xt(tz) = 0.85. Solving the minimization equations (28) to find tan/3 
and /it, will then determine all remaining parameters of the model. Let us first try to assess 
the accuracy of the approximation X t ^> A&, A r that we shall use. This is expected to 
be a reliable one provided tan/3 is of order unity. Thus, starting, say, with tan/3 = 5.5 

and Y x (tz) = [5.810 -3 , 2.010 -4 , 6.410 -5 ], [x = u,d,e] and integrating upwards to large 
scales in the above approximation, yields: Y x (0) = [3.710 -4 , 1.210 -5 , 7.710 -6 ]. Inte- 
grating now downwards from Mx with the exact evolution equations, yields: Y x (tz) = 
[3. 110 -3 , 1.610 -4 , 2.210 -5 ], which is in fair agreement with the starting values. 

The numerical applications presented in this work should usefully complement those 
of ref.[19]. The Calabi-Yau and O-II (orbifolds with small threshold corrections) scenar- 
ios have no impact on the flavor changing neutral current issue. Building up on the 
(large threshold correction) O-I scenario of ref.[19], we consider here a wide variety of non- 
universal modular weights and also include the b-quark and r-lepton Yukawa couplings. 
The latter item opens up the possibility of independently fixing tan/3 when solving for 
the radiative electroweak breaking condition. However, the approximation At ^> Aft,A r , 
eq.(19), that we use is a limiting factor here, because this weakens the sensitivity to these 
coupling constants of the electroweak symmetry breaking equations. We should mention at 
this point that our method of solving for tan/3 is not very robust. We proceed numerically 
along the following steps: First, we select a trial value for tan/3, which fixes Xb(tz), X T (tz)- 
Second, we evaluate all other parameters by solving the renormalization group equations 

for a discrete set of values for ^^- inside the interval [1 — 4]. Finally, we solve for /i(0) 

the combination of the two equations (28) obtained by elimination of the explicit depen- 
dence on tan/3; we determine the associated value of tan/3 and run down again the scale 
evolution equations with these up-dated values of //(0) and tan/3. While this procedure 
leads to accurate solutions for the second equation (28), thanks to its strong sensitivity to 
/i(0), it fails in general to accurately solve the first equation. The inaccuracies originate 
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from the disparity between the two scales in the problem, m g and rriz, and worsen with 
increasing m g . For m g < 200 GeV, they represent from a few 10 % in favorable cases to 
50 %. However, the inaccuracies dramatically increase with larger values of m g . 

Results are shown in figure 1 for case I-a and figure 2 for case Ill-a. We see (figures a-b) 
that solutions with tan/? of order unity and concomitantly //(0) large in comparison to m g , 
are generally selected. The variation with 9 of these solutions to the electroweak breaking 
equations depend on the boundary conditions used for Higgs bosons masses. For widely 
unequal Higgs bosons masses at Mx (cases I or II), with increasing 9, tan/3 decreases and 
\x increases; the opposite takes place for equal masses at Mx (case III). The physical (scale 
Tfiz) masses of gauginos, scalars and fermionic superpartners (figures (c), (d), (f)) are all 
increasing functions of | sin#|. For gaugino masses, this is due to the fact that the dilaton 
term dominates over the moduli term. While the latter term increases with < T > and 
with 5gs, in all cases that we deal with here, it represents a small, few % fraction of the 
former and is of negative sign. Note that the moduli contribution is what prevents mass 
ratios of the three gaugino masses from being ^-independent. 

The scalars mass matrices in the interaction basis are diagonal ones, due to our re- 
striction to canonical kinetic terms and the simplified version used for the scale evolution. 
The values of the averaged traces of the scalars mass matrices are seen to lie at a factor 
2 — 4 (1) times the basic energy scale, m g for squarks (sleptons). The third generation 
masses only are affected by the Yukawa couplings. These contributions shift down the 
(3, 3) scalars masses in comparison to the (1, 1), (2, 2) masses by a few % for squarks and 
10 — 20% for sleptons. (Since we do not exactly diagonalize the full mass matrices, by 
taking into account left-right chirality mixing, the foregoing statements are approximate 
ones.) The top-quark and Higgs boson masses, which have not been plotted in the fig- 
ures, are also smooth increasing functions of | sin#|. As 9 sweeps the interval [0, n], the 
top-quark mass in case III, for instance, increases monotonically in the interval [111 — 150] 
GeV; the corresponding interval for the tree level Higgs boson mass is [12 — 90] GeV. For 
case I, the variation with 9 is of opposite sense. The third generation trilinear parameters 
A x (tz) (figure (e)) are strongly influenced by radiative corrections. For case I, as 9 sweeps 
the interval [0, n], A t (0) regularly decreases from 1.4 to —3.4, while A t {mz) (figure 1 (e)) 
first levels off then slowly decreases from 4 to 0.6. The parameters Ab )T are less affected by 
radiative corrections. This is also the case for the B— parameter, for which -B(O) regularly 
decreases from —0.4 to —3.6 while B(tz) decreases from —0.3 to —2.5. 

The rapid decrease with decreasing 9 of the average scalars masses near small 9 or n — 9 
is explained by the fact that the masses of scalars with finite modular weights \n a \ > 1, 
may then take negative values at Mx- When these values are large enough in absolute 
value, as for the case of large \n a \, the repulsive contributions of gaugino interactions, 
which also get smaller there because of the reduced gauginos masses, are unable to flip 
the physical scalars masses sign to positive. This results in drastically reduced averaged 
traces. Such values of 9 should be discarded since the associated vacuum would then 
lead to tachyonic particles. The main agent to avoid such negative squared masses are 
gaugino loops. These are more effective for squarks than for leptons. The contributions 
of Yukawa couplings oppose those of gauginos and are obviously larger for u-squarks than 
d-squarks and sleptons. As a result, the upper bounds on cos 9 from the constraint of 
the absence of tachyonic particles are mainly determined by the sleptons masses. These 
bounds depend somewhat on the modular weights assignments. For — or |— — 1| > 0.2 no 
tachyonic particles are present in nearly all three cases described above. This is the reason 
why we have deleted these intervals in plotting the figures. Finally, we observe that the 
chargino (neutralino) minimal mass eigenvalues (figure (f)) lie at 3 (|) times m g with a 
#— dependence similar to that of the wino mass. 
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The parameter m g essentially fixes the overall scale for most dimensionful quantities 
in spite of the fact that the ratio — - is also relevant to the electroweak breaking condition. 

Increasing m g in the interval [100 — 1200] GeV, affects slightly tan/3 and negligibly the 
dimensionless parameters, A X ,B, while all superpartner mass parameters essentially scale 
linearly with m g . This is clearly demonstrated by figure 3. 

Let us briefly describe the effect of changing the interval of variation of the angle 
9 from the first and second quadrants, corresponding to the choice made above, to the 
third and fourth quadrants. This change leaves scalar masses unchanged but modifies the 
interference of moduli and dilaton contributions in the other parameters, flipping the signs 
of A x , B and leading to gaugino masses of negative sign. Naturally we need to flip the 
signs of gaugino masses to positive by performing a phase redefinition of the Majorana 
spinor fields of gauginos. The results for the interval of angles 9 > n are not qualitatively 
different from those above in the interval 9 < n. 

3.5 Flavor changing parameters 

The flavor changing parameters are dimensionless and hence essentially independent 
of m g . Their variation with Goldstino angle is depicted in figure 1 for case I-a and figure 2 
for case Ill-a (figures (g)-(l)). Let us observe first that because absolute values are plotted 
here, the sudden changes of slopes featured by certain curves are simply due to changes of 
sign in the corresponding flavor changing parameters. The frequent occurrence of such sign 
flips, as a function of 0, reflects the presence of large cancellations between the boundary 
conditions and the radiative corrections contributions which vary in opposite sense with 
variable 9. 

Considering first the unmixed chirality parameters, we see that they are decreasing 
functions of |sin#|; sleptons have the fastest slopes and d-squarks come next. The fast 
decrease near small 9 or n — 9 is a combined effect of the decreasing averaged scalars 
masses there and the increasing generation dependent boundary conditions (cf. eqs.(36), 
(40)). Due to the specific structure of the fermions transformation matrices, the parameters 
(Smm)ij are roughly proportional to the differences of modular weights rii — rij. Turning 
next to the mixed chirality parameters, we see that they vary rather slowly with 9. This is 
explained by the fact that the non-universal contributions vary in an opposite sense from 
the gaugino contributions, as previously discussed in subsection 3.2. Sleptons have the 
largest slopes because of the smaller gauginos contributions there. 

A comparison with experimental bounds in Table 1 reveals that our predicted values 
in figures 1 and 2 lie safely below the experimental ones, with the exception of 8f 2 an d 
(^lr) 12 f° r x = d,e. We can infer from this comparison lower bounds on 9 and n — 9. For 
case I-a, the lower bounds in the d-squarks case lie approximately at 0.2, while the entire 
interval of 9 angles is excluded in the sleptons case. For case Ill-a, the squarks case entails 
a larger lower bound of 0.4, while the sleptons case again excludes all angles. The excluded 
intervals in 9 increase in proportion to the size of perturbations away from universality 
of modular weights and also to the size of modular weights themselves. For the mixed 
chirality parameters, the model predictions largely exceed experimental bounds, except for 
(^lr) 12 f° r which predictions lie an order of magnitude above the bounds for the entire 9 
interval. 

For a refined comparison one needs to rescale experimental bounds in Table 1 by 
factors (m/lOOOGeV) for squarks and (m/70GeV) p for sleptons, with p = 2 for 8%j M 
and p = 1 otherwise, using the predicted average scalar masses given by figures (d). 
This correction induces on the squarks bounds a reduction by factors of average size 2 
and for sleptons an enhancement by factors of average size 10 (3) for unmixed (mixed) 
chirality parameters. This will reduce the predictions for squarks, hence will increase the 
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lower bounds on 9 and n — 9, but will partially fill in the gap between predictions and 
experimental bounds for sleptons. The rescaling with respect to the variable x is not very 
important for squarks, for which the model predictions are compatible with x ~ 1. For 
sleptons, however, a wide gap separates the values occurring in predictions, x 1 ~ 1, from 
that used in Table 1, x 1 ~ 0.02. Since a rescaling of the experimental bounds in Table 1 
to x 1 = 1 is expected to raise these bounds by a factor of order unity [34] , this would not 
modify the above conclusions. We have not pursued here a more quantitative study of the 
x 1 dependence. 

It is interesting to determine the maximal perturbations away from universality of 
modular weights allowed by the condition of absence of tachyonic particles and the bounds 
on flavor changing parameters in the first and second generations. We find that perturbing 
the solutions in ref.[8] for modular weights by a generation non-universality exceeding 
two units results in flavor changing parameters larger than experimental over the entire 
interval of 9. This would essentially rule out the model. Detailed results for various cases 
are presented in Table 2. The three versions of case I are characterized by degenerate first 
and second generations. The one unit differences of modular weights with respect to the 
third generation results in contributions of size 10 -3 , hence roughly the same order as the 
radiative corrections (compare I-a and I-a'). The comparison of results between cases I-a 
and I-a' and Il-b and II-b' shows us that the radiative corrections from first and second 
generations may contribute at about the same level as the boundary conditions. Certain 
changes of signs reflect the presence of strong cancellations between radiative corrections 
and boundary conditions contributions. The results for case II show that mixed chirality 
are proportional to difference of modular weights , since going from Il-a to II-b enhances 
the parameters by factor 4 — 10. The right chirality parameters are more strongly enhanced 
because of stronger cancellations effects present there. Finally, we note that the mixed 
chirality parameters are remarkably stable within factors 2 — 3. 

In order to reduce the flavor changing parameters two main possibilities can be en- 
visaged. One can eventually assign untwisted string modes, with n a = —1, to all three 
generations or to a subset. This choice might conflict with the requirement of obtaining 
large scale hierarchies in the fermionic mass matrices. An alternative possibility is to allow 
for a non-vanishing goldstino matter component. This could significantly reduce the non- 
universality, since, as seen from eq.(10) the non universal contributions involve a cos 2 6 a 
factor reduction which is stronger than the opposing gauginos loops universality dilution 
involving cos 9 a- In this case one would reduce both squarks and sleptons parameters. 

4. Conclusions and outlook 

The main issue adressed in this paper concerned the flavor changing neutral current 
constraints imposed on standard-like orbifold models using the duality symmetry approach 
initiated in refs.[8,19]. We have obtained detailed numerical predictions for the low energy 
spectrum and flavor changing parameters, essentially covering the entire parameter space 
of the model. Our results complement those of ref.[19] and establish that all experimental 
bounds for masses of superpartners scalars and fermions can be obeyed provided the grav- 
itino mass is chosen to lie above roughly 200 GeV. While the flavor changing parameters are 
controlled at unification scale by a strong dependence on modular weights and goldstino 
angle the residual contributions left out after radiative corrections are taken into account 
are proportional to modular weights (differences) and have a smoothed dependence on 
9. The relevant flavor changing parameters in the comparison with existing experimental 
bounds are those associated with the first and second generations, especially sleptons pa- 
rameters. The experimental bounds can be comfortably respected if 9 or n — 9 lie above 
0.3 and if non-universal perturbations are chosen for modular weights not exceeding two 
units away from values in solutions compatible with anomalies cancellation and gauge cou- 
pling constants unification. Lower bounds of same size are also required in order to avoid 
tachyonic particles. The main features of the model which imply these properties are the 
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large gauginos and scalars masses in comparison to the gravitino mass. The framework 
in which the present work has been developed involves three main shortcomings which we 
briefly review here in order to motivate a future improved treatment of the problem. 

(i) Off-diagonal couplings in the matter Kahler potential. Terms of form, SK = 

f a p(M : M)A a A*g + c.c, with a, f3 representing identical gauge group low energy modes but 

in different generations, are not excluded in general, although these are generically absent 
in orbifold models. Even when present, such terms can always be transformed away by a 
linear transformation of superfields, which however will affect the Yukawa couplings. Off- 
diagonal matter Kahler potentials can also arise if massless string modes mix, so that low 
energy fields are identified with linear combinations of string modes. Such a situation is 
commonly encountered when one exploits flat directions to assign VEVs to certain scalars 
designed to reduce the rank of the gauge group or to truncate the massless spectrum. 

(ii) Relative alignment of scalars mass and Yukawa coupling matrices in generation 
space. On side of gaugino loops radiative corrections, this is the next important item 
needed to explain naturally small flavor changing parameters. Although the structure of 
the fermions and scalars superpartners matrices are sensitive to different characteristics 
of orbifolds ( modular weights versus fixed points) it would be desirable to demonstrate 
on specific examples whether some relationship between them exists and to demonstrate 
its properties. While the threat from a minimal non-universality which is restricted to 
diagonal elements of mass matrices has been partially discarded through our present study, 
one still needs to worry about off-diagonal terms. Such contributions are difficult to avoid 
when non-renormalizable interactions are taken into account. It is then necessary to verify 
whether the occurrence of large hierarchies in scalars mass matrices and in Yukawa coupling 
matrices are compatible and if this can be achieved in orbifold models without postulating 
an horizontal symmetry. 

(Hi) CP-violation bounds. Well-known mechanisms contributing to imaginary parts of 
mass shifts in the neutral mesons systems or to electric dipole moments provide estimates 
for the imaginary parts of flavor changing parameters which are an order of magnitude 
smaller than the experimental bounds for real parts. Approximate calculations of the 
relevant CP-violating phases, <pA = Im(AM a ), 4>b = Im(Bff), indicate the need of 
assigning small values to the complex phases in the scalar or auxiliary components of 
superfields. If such small contributions are confirmed by a more quantitative study, this 
would then require that dilaton and moduli fields VEVs should stabilize very close to 
real values, a property which seems to be realized in a natural way in descriptions of 
supersymmetry breaking incorporating duality symmetries. 
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Tables Captions 

Table 1. Experimental upper bounds for the S MN matrices of up-squarks (1st line), d- 
squarks (2d-4th lines) and sleptons (5th line). Our sources are the works of Gabbiani 
and Masiero (GM) [20], Nir and Seiberg (NS) [27] and Hagelin et al., (HKT) [341. We 
set the average squarks mass to m = 1000 GeV and sleptons mass to to = 70 GeV. The 

parameters x 3 = (^ff) 2 (squarks) or x 1 = (~^-) 2 (sleptons) are set to 1 in the bounds of 
NS [27] and HKT [34] and to x 3 = 0.49, x 1 w 0.029 in those of GM [20]. (Contributions 
arising from A and \x terms in ref.[20] are neglected.) The notation 5ij is reserved to bounds 
obtained from box diagrams contributions proportional to the products SllSrr- Likewise, 
bounds on 5lr are obtained from contributions involving the products [(5 LR)ij{.& lr) j^ 11 • 
For neutral KK, BB.DD systems, the bounds on (Smn)^ scale as to at fixed x. For 
decays qi — ► qj + 7, li — ► lj + 7, the bounds for (Smm)ij scale as to 2 and those for (Slr)^ 
as to, at fixed x. Thus, except for the sleptons bounds 5mm here, which scale as to 2 , all 
other bounds quoted in the table scale as to. Note that the new Cleo measurements for 
b -> s + 7 yield the bound (<J£ fl ) 2 3 = 0.028. 

Table 2. Flavor changing parameters for d-squarks (top half) and sleptons (bottom half) 
for fixed m g = 2QQGeV, 9 /n = 0.35. The prime indicates that we have omitted radiative 
corrections in first and second generations which are described in the leading logarithm 
approximation of eq.(23). The various cases described in Section 3.3 are identified in 
the first column. The modular weights in the order a = [Q; U c ; D c ; L; E c ] are: case 
I-a: n a = -[221; 332; 221; 223; 223]; case Il-a: n a = -[111; 133; 111; 333; 133]; casell-b: 
n a = -[023; 133; 023; 023; 133]; case III-a:n a = -[123]; case III-b:n Q = -[112]; 



Figures Captions 

Figure 1. Plot as a function of 6/ V for m g = 200 GeV of the model parameters 
in Case I-a: (a) tan/?; (b) /u(0); (c) standard model gauginos masses M a (tz) in the 
order a = 3, 2, 1 from top to bottom ; (d) Average trace of left times right chiral- 

ity mass submatrices, fh x = [Trace(M'[ L (tz))Trace(M] lR (tz))] 1 , for d-squarks (con- 
tinuous), u-squarks (dot-dashed) and sleptons (dashed); (e) supersymmetry breaking 
parameters A x (tz) for x=t (continuous), b (dot-dashed), r (dashed), and B(tz) (dot- 
ted); (f) minimum masses of charginos (continuous) and neutralinos (dot-dashed); (g)-(i) 

S fj = \( S LL)ij( S RR)ij)\^ \ x = d,e,u] for (i,j) = (1,2) (figure (g)), (1,3) (figure (h)) and 
(2,3) (figure (i)); (j)-(l) % = W lr)ij {51 r ) j% )\^ [x = d,e,u] for (i,j) = (1,2) (figure 
(j)), (1,3) (figure (k)) and (2,3) (figure (1)). Curves in each of figures (g)-(l) refer to d- 
squarks (continuous), sleptons (dot-dashed) and u-squarks (dashed). Note that in figures 
(g)-(l) the discontinuities in slopes occurring for certain curves signal the occurrence in 
these semi-logarithmic plots of changes of signs for quantities in ordinate. 

Figure 2. Same conventions as figure 1 in Case Ill-a. 

Figure 3. Plot versus m g for fixed 6 = 0.357T of low energy parameters in Case I-a, using 
the same conventions for curves as in figure 1. 
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TABLE 1 



Sector (5 m m) 12 (5 mm) 13 (5 mm)23 812 ^13 #23 (5 lr) 12 ($lr) 13 (£1^)23 

U (NS) 0.1 0T04 (L06 

ref.[27] 

D (NS) 0.05 01 0.006 O04 0.008 (L06 51)4 

ref.[27] 

D (GM) 0.035 0.23 0.0051 0.074 0.007 0.1 0.12 

ref.[20] 

D (HKTp.l 0.27 47 0.006 0.073 0.0044 0.14 0.071 

ref.[34] 

L (GM) 8.9610" 4 7.67 708 1.5510- 6 2.2410- 1 2.0610" 2 

ref.[20] 
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TABLE 2 



Case 


($ll) 12 


(Sll)iz 


(Sll)23 


(Srr) 12 


(Srr)i3 


(Srr) 23 


(Slr)i 


12 


(Slr): 


13 


(^Li?)23 


D-SQUARKS 

la 3.210-3 
la' -4.010" 
Ib 3.210-3 


-3.610" 
4 -5.110" 
-8.510" 


-13.010" 3 

-13.710-3 

t^lO" 3 


3.710" 4 
-4.310" 
3.310" 5 


4.710" 4 
t^lO" 4 
4.110" 4 


-6.610" 
-5.110" 
-4.210" 


19.210" 
19.210" 

t.gio- 


-5 
-5 
-5 


9.410" 
9.410" 
8.410" 


-5 
-5 
-5 


4.210" 4 
4.210" 4 
3.710" 4 


Ha 

lib 
lib' 


2.910-3 
1.110" 2 
7.310-3 


-6.410" 

-9.710" 
-1.110" 


-^lo- 3 

" 4 -5.810" 
-3O.610-3 


6.710" 5 
38.410" 3 
7.810" 3 


3.110" 4 
1.310" 4 
1.110" 4 


6.810" 5 
8.810" 4 
9.710" 4 


9.810" 
1.210" 
1.210" 


-5 
-4 
-4 


1.110" 
1.210" 
1.210" 


-4 
-4 
-4 


4.710" 4 
5.210" 4 
5.210" 4 


Ilia 
Illb 


5.810-3 
1.910-3 


-1.310" 
-1.210" 


"37.210-3 

-3o.810-3 


3.010-3 
-2.310" 


1.610" 4 
t.610" 4 


7.210" 4 
7.210" 4 


9.210" 
8.510" 


-5 
-5 


1.010" 
1.010" 


-4 
-4 


4.510" 4 
4.610" 4 


SLEPTONS 

la 7.110-3 
la' 2.310" 4 
lb 1.010" 2 


2.810" 4 
-1.510" 
4.110" 4 


-7.710" 

-t.sio" 3 

-2.810" 


4 -1.710" 
-2.210" 
3-1.810" 


- 4 -4.810" 
- 4 -4.810" 
- 4 -6.910" 


l.llO" 2 
l.llO" 2 
\A10~ 2 


3.410" 
3.410" 
3.610" 


-5 
-5 
-5 


4.110" 
4.110" 
4.310" 


-5 
-5 

-5 


6.010" 4 
6.010" 4 
6.310" 4 


Ha 
lib 
lib' 


7.010-3 
2.210" 2 
1.410" 2 


6.810" 4 
3.410" 4 
-1.310" 


-6.110" 
-1.810" 
t.OlO" 3 


32.610" 2 
^lO" 2 
2.610" 2 


3.110" 4 
3.110" 4 
3.010" 4 


3.710" 5 
6.110" 5 
1.110" 4 


4.710" 
5.510" 
5.510" 


-5 
-5 

-5 


3.210" 
3.310" 
3.310" 


-5 
-5 
-5 


4.510" 4 
4.910" 4 
4.910" 4 


Ilia 
Illb 


1.810" 2 
6.910-3 


3.510" 4 
2.110" 4 


4.710" 4 
2.210" 3 


1.810" 2 
-2.010" 


-7.510" 
- 4 -6.410" 


l.sio- 2 

I.SIO" 2 


4.810" 
3.110" 


-5 
-5 


4.010" 
3.510" 


-5 
-5 


6.010" 4 
5.410" 4 
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